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Abstract 

The most important pathogenic fungus of the silkworm, Bombyx mori L. (Lepidoptera: Bombyci- 
dae), is Beauveria bassiana (Balsamo-CrivelU ) Vuillemin (Hypocreales: Clavicipitaceae), which 
causes significant damage to sericulture production. Therefore, diagnosing fungal disease and 
developing new control measures are crucial to silk production. To better understand the respon- 
sive and interactive mechanisms between the host silkworm and this fungus, variations in 
silkworm gene expression were investigated using the suppression subtractive hybridization 
method following the injection of B. bassiana conidia. Two cDNA libraries were constructed, 
and 140 cDNA clones were isolated. Of the 50 differentially expressed genes identified, 45 (112 
clones) were identified in the forward library, and 5 (28 clones) were identified in the reverse li- 
brary. Expression profiling of six of these genes by quantitative polymerase chain reaction 
(qPCR) verified that they were induced by the fungal challenge. The present study provides in- 
sight into the interaction between lepidopteran insects and pathogenic fungi. 
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Introduction 

Beauveria bassiana Balsamo-Crivelli Vuille- 
min (Hypocreales: Clavicipitaceae) is an 
entomopathogenic fungus that grows naturally 
in soil and has been widely used in the biolog- 
ical control of insects in forests. However, its 
mechanism of action remains elusive despite 
chemical and synthetic methods of analysis 
(Grogan and Holland 2000; Hao et al. 2001). 
In the natural environment, B. bassiana usual- 
ly infects insects and causes white muscardine 
disease through cuticle penetration. The fun- 
gus is then fought by the innate immune 
responses of the insects, including cellular and 
humoral mechanisms (Lavine and Strand 
2002; Hou et al. 2011). Entomopathogenic 
fungi replicate in the insect hemolymph. 
However, the insect genes that respond to 
fungal infection and replication are still not 
well known. 

The silkworm, Bombyx mori L. (Lepidoptera: 
Bombycidae), is a typical lepidopteran insect 
and is economically important for silk produc- 
tion in many developing countries. It has also 
contributed enormously to the study of insect 
genetics and immunology (Goldsmith et al. 
2005; Li et al. 2005; Liu et al. 2009). 

B. bassiana is one of the major fungal patho- 
gens of the silkworm and can cause enormous 
damage to the sericulture industry. The mo- 
lecular mechanism by which it infects 
silkworms is still poorly understood. There- 
fore, elucidating the mechanism of antifungal 
immunity of the silkworm is important to im- 
prove its antifungal ability. The identification 
and study of differentially expressed insect 
genes after immune challenge will lead to a 
better understanding of the breadth and regu- 
lation of insect immune responses. 
Suppression subtractive hybridization (Di- 
atchenko et al. 1996) has been used to identify 
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immune-inducible genes in the mosquito 
(Oduol et al. 2000), tsetse fly (Hao et al. 
2001), and Manduca sexta (Zhu et al. 2003). 
In the present study, suppression subtractive 
hybridization was used to identify the differ- 
entially expressed genes between Beauveria- 
injected and water-injected silkworm larvae, 
and these data were verified by real-time 
quantitative polymerase chain reaction 
(qPCR). 

Materials and Methods 

Silkworm strain 

The silkworm strain Dazao, provided by the 
Sericultural Research Institute of the Chinese 
Academy of Agricultural Sciences, was used 
in the study. The larvae were reared on fresh 
mulberry leaves at 25° C. Third-day fifth 
instar larvae were used in the experiments. 

Treatment with B. bassiana conidia 

B. bassiana conidia were diluted to a 

o 

concentration of 10 spores/mL with sterile 
distilled water. Each larva was injected with 1 
|iL of the conidia solution. The larvae of the 
control group were injected with the same 
amount of sterile distilled water. The larvae 
were reared under high temperature and high 
humidity conditions (28° C and 95% RH) for 
8 hours to promote conidia germination and 
then returned to normal conditions (25 °C and 
80% RH). 

Cuticle and hemolymph collection 

After injection with the conidia and sterile 
distilled water, the cuticles and hemolymph of 
the larvae were collected at 3, 6, 9, 12, and 24 
hours after injection. The hemolymph was 
directly mixed with pre-joined Trizol reagent 
(Life Technologies, 
www. lif etechnologies . com ) in an Eppendorf 
tube. The cuticles were quickly washed in a 
diethylpyrocarbonate-treated phosphate 
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buffered saline solution (137 mM NaCl, 2.68 
mM KCl, 8.1 mM Na2HP04, and 1.47 mM 
KH2PO4 at pH 7.4) and immediately frozen in 
liquid nitrogen. 

Isolation of total RNA and polyA+ RNA 

The total RNA from the cuticles and 
hemolymph of the B. bassiana-'m]QciQd and 
water-treated larvae was extracted using 
TRIzol reagent and subjected to DNase I 
(Takara Bio, www.takara-bio.com ) treatment 
according to the manufacturer's instructions. 
The total RNA concentration was determined 
using a Biophotometer (Eppendorf, 
www.eppendorf.com ) by measuring the 
absorbance values at 260 nm and 280 nm 
(A26o-A28o). Equal amounts of total RNA from 
the samples of the five time points were 
pooled. The Oligotex mRNA Mid Kit 
(Qiagen, www.qiagen.com) was used to 
purify polyA+ RNA from the RNA pool, and 
2 \ig of poly(A)+ RNA was used as the 
starting material for reverse transcription to 
construct the subtracted cDNA libraries. 

Construction of cDNA libraries by 
suppression subtractive hybridization 

Suppression subtractive hybridization was 
performed using a Clontech PCR-select 
cDNA subtraction kit (Clontech Laboratories, 
www. clontech. com ) . Reciprocal forward and 
reverse subtractions were performed 
according to the manufacturer's instructions. 
The RNA extracted from the cuticles and 
hemolymph at all the time points was equally 
mixed. The forward library was constructed 
using the cDNA of the B. bassiana-'m]QciQd 
cuticles and hemolymph as the tester and the 
cDNA of the control as the driver. The reverse 
library was constructed using the cDNA of the 
control as the tester and the cDNA of the B. 
bassiana-'m]QciQd cuticles and hemolymph as 
the driver. 



Hou et al. 

The subtracted cDNA libraries were generated 
by inserting the differentially expressed 
cDNA fragments into the pGEM-T Easy vec- 
tor (Promega, www.promega.com ) and 
transforming these vectors into JM109 compe- 
tent cells. Aliquots (100 |iL) of the 
transformation mixture were then spread onto 
Luria-Bertani agar plates containing 100 
mg/mL ampicillin, 80 mg/ml X-gal, and 50 
mM isopropyl 1 -thio-(3-D-galactopyranoside 
and incubated at 37° C overnight. Some sub- 
tractive clones were sequenced, and the 
nucleotide and amino acid sequence homolo- 
gies were determined by searching the 
NCBI/GenBank database using the BLASTX 
algorithm. 

Confirmation of differentially expressed 
genes by real-time qPCR 

The total RNA from the cuticles and 
hemolymph of B. bassiana-infQctQd and 
control larvae was extracted using Trizol 
reagent (Life Technologies). The RNA was 
treated with DNase I following the 
manufacturer's instructions. The 

concentration of DNase I-treated RNA was 
adjusted with H2O to 1 |ig/|iL, and 1 |ig of the 
DNase I-treated RNA was reverse transcribed 
in a 20 |iL reaction using a PrimeScript RT 
reagent kit (TaKaRa). Real-time qPCR was 
performed using 2 |iL of the diluted first- 
strand cDNA (1/100) in each 20 |iL reaction 
mixture. The reaction was performed with 
specific primers for amplifying the following 
genes: chemosensory protein 11, muscle LIM 
protein isoform 1, transferrin, arylphorin, sex- 
specific storage-protein SPl precursor, 
lysozyme, moricin, and the low molecular 
lipoprotein 30K pBmHPC-6 (Lp-c6) (Table 
1). The relative expression levels of each gene 
at different time points were normalized using 
the Ct values obtained for B. mori 
glyceraldehyde-3 -phosphate dehydrogenase 
(5mGAPDH), an endogenous control gene. 



Journal of Insect Science | www.insectscience.org 



3 



Journal of Insect Science: Vol. 13 | Article 138 

amplification in the same plate. In each assay, 
the expression level is shown relative to the 
lowest expression level, which was arbitrarily 
set to one. All samples were tested in 
triplicate. The mean values and the standard 
deviations were used for the analysis of the 
relative transcript levels for each time point 
using the relative quantitative method (AACt). 

Real-time qPCR was performed using 2 |iL of 
diluted first strand cDNA (1/100) in each 20 
|iL reaction volume using SYBR Premix Ex 
Taq (TaKaRa) according to the manufactur- 
er's instructions. Specific primers for the 
genes chemosensory protein 11, muscle LIM 
protein isoform 1, transferrin, arylphorin, sex- 
specific storage-protein SPl precursor, low 
molecular lipoprotein 3 OK pBmHPC-6 (Lp- 
c6), lysozyme, moricin, and 5mGAPDH are 
listed in Table 1. The final concentration of 
each primer was 100 nM. The reactions were 
run in triplicate on an Opticon system (Bio- 
Rad, http://www.bio-rad.com ) using the 
following thermal cycling parameters: 95° C 
for 5 sec and 40 cycles at 60° C for 10 sec and 
72° C for 10 sec. The melting curves were 
constructed after amplification. The data were 
analyzed and normalized relative to the 
5mGAPDH transcription levels using Opticon 
monitor analysis software (BioRad). The 
AACt method was used to evaluate the relative 
expression differences. 

Results 

EST sequencing and identification 

The cDNA libraries, both forward and re- 
verse, were constructed from the silkworm 
larvae injected with B. bassiana conidia using 
suppression subtractive hybridization. The 
forward library was constructed with cDNA 
from the B. bassiana-'m]QciQd cuticles and 
hemolymph as the tester and cDNA of the 
control as the driver, whereas the reverse li- 
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brary used cDNA of the control as the tester 
and cDNA of the B. bassiana-'m]QciQd cuticles 
and hemolymph as the driver. A total of 140 
cDNA clones were isolated, 112 from the 
forward library and 28 from the reverse li- 
brary. In total, 45 genes were identified in the 
forward library, and five genes were identified 
in the reverse library (Tables 2 and 3). 

The genes present in the forward library were 
upregulated by B. bassiana infection (Table 2) 
and were more abundant than the genes in the 
reverse library (Table 3). Using blast2G0 
software and the annotations of the Spodop- 
tera frugiperda sequences (Barat-Houari et al. 
2006), the subtractive genes were classified 
into five groups: genes encoding proteins 
ubiquitously expressed by many cell types 
(AI-AIX), genes responsible for cell-cell 
communication (BI-BIII), genes encoding 
transcription factors and gene-regulatory pro- 
teins (C), genes encoding molecules expressed 
in insects (DI-DIV), and others (EI-EIII). 
There was only one overlapping gene between 
the two libraries, which encodes the B. mori 
30 kDa lipoprotein 19G1 precursor (Tables 2 
and 3). 

Some upregulated genes identified in our pre- 
vious study on the percutaneously infected 

silkworm were also induced in the present ex- 
periment, such as putative cuticle proteins, 
heat shock proteins, ribosomal proteins, and 
the antimicrobial proteins lysozyme and the 
cecropin B precursor. The newly identified 
genes in the present study included those en- 
coding ferritin, transferrin, ATP synthase 
subunit, troponin, storage proteins (arylphorin, 
sex-specific storage protein (SPl), and sex- 
specific storage protein 2 (SP2)), a 30 kDa 
lipoprotein, and chemosensory protein 11 
(CSPl 1). The genes that encode ferritin, trans- 
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ferrin, heat shock proteins, CSPll, and stor- 
age proteins were immune-responsive, and the 
genes that encode moricin and gloverin 4 were 
anti-microbial, whereas the genes that encode 
ribosomal proteins, elongation factors, and 
ATP synthase were involved in the pathologic 
processes. Furthermore, ribosomal proteins 
also participated in signal transduction. 

Determination of differentially expressed 
cDNAs by real-time qPCR 

The differential gene expression in the larval 
cuticles and hemolymph upon B. bassiana in- 
fection was determined via real-time qPCR 
(Figure 1). The present study focused on the 
genes that were present in the forward library, 
as they are likely responsive to the B. bassi- 
ana infection and may contribute to the 
resistance of the silkworm against B. bassiana 
infection and propagation. The transcript 
abundances of the selected eight genes ex- 
pressed differently in the hemolymph between 
the B. bassiana-'m]QciQd and the water- 
injected larvae were compared using real-time 
qPCR (Figure 1). 

In the hemolymph of the B. bassiana-'m]QciQd 
larvae, the transcript level of the CSPll gene 
was the highest at 6 hours after injection and 
then gradually decreased, whereas the tran- 
scription in the control larvae stayed almost 
constant in the same period (Figure lA). The 
hemolymph transcript levels of the low mo- 
lecular lipoprotein 3 OK pBmHPC-6 (Lp-c6) 
gene exhibited the same expression pattern for 
both B. bassiana-'m]QciQd and water- injected 
larvae and reached a peak at 6 hours after in- 
jection; however, its transcript levels in the B. 
bassiana-'m]QciQd larvae were much higher 
than those in the water-injected larvae (Figure 
IF). The gene transcript levels of muscle LIM 
protein isoform 1, transferrin, sex-specific 
storage-protein SPl precursor, and arylphorin 
reached their peaks at 9 hours after injection 
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and then gradually decreased in both the B. 
bassiana-'m]QciQd and the water-injected lar- 
vae. However, the transcript levels in the B. 
bassiana-'m]QciQd larvae were much higher 
than those in the water-treated larvae, and the 
transcript level of the arylphorin gene was the 
highest among the six genes (Figure IB-IE). 
In the B. bassiana-'m]QciQd larvae, the tran- 
script levels of lysozyme and moricin reached 
their peaks at 6 hours and 9 hours after injec- 
tion, respectively, and then gradually 
decreased; however, their expression levels 
not significantly different from the water- 
injected control. The transcription of these 
genes in the cuticles in the same period was 
also detected, but the differential expression 
between the B. bassiana-'m]QciQd and water- 
injected larvae was not obvious and almost 
negligible (data not shown). 

Discussion 

The mechanisms of immune responses to bac- 
teria in silkworms have been extensively 
studied. The responsive mechanisms of in- 
sects to fungal infection and some antifungal 
peptides have drawn the attention of research- 
ers, and related signal transduction molecules 
have been identified (Belvin and Anderson 
1996; Hoffmann 2003). Differential gene ex- 
pression in silkworm larvae percutaneously 
infected with B. bassiana was reported in our 
previous study (Hou et al. 2011). To better 
understand the silkworm responsive mecha- 
nism to fungal infection, silkworms were 
infected in the present study by injecting B. 
bassiana conidia into the hemolymph of the 
larvae, and the genes that were differentially 
expressed upon infection were analyzed. 

After fungal conidia enter the hemolymph of 
the silkworm, a series of immune changes oc- 
curs in the hemolymph (Hoffmann and 
Reichhart 2002; Tzou et al. 2002). The ex- 
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pression levels of some genes that are respon- 
sive to fungal invasion also change after 
fungal injection, as identified in the present 
study. Most of these genes are functional; 
their expression was upregulated, and they are 
likely to play important roles as the infection 
progresses. In the present experiment, only a 
few differentially expressed genes were de- 
tected in the reverse library. These results are 
similar to those of our previous study (Hou et 
al. 2011). Therefore, the subsequent analysis 
focused on the identification of differentially 
expressed genes present in the forward library 
because these genes are likely responsive to B. 
bassiana infection and may affect B. bassiana 
proliferation in the silkworm larvae. 

The differentially expressed genes identified 
in the forward library were upregulated, which 
implies that these genes are induced by or re- 
sponsive to B. bassiana infection. Some of 
these genes encode storage proteins, such as 
arylphorin, SPl, SP2, and some 30 kDa lipo- 
proteins. These storage proteins are the major 
nutrient sources for growth and development, 
especially for insect metamorphosis (Willott 
et al. 1989; Spyliotopoulos et al. 2007). They 
are also constituents of the sclerotizing system 
of the cuticle (Dong et al. 1996). Among these 
storage proteins, arylphorin transcript levels 
were suppressed during parasitization of Heli- 
othis virescens larvae by Campoletis 
sonorensis (Shelby and Webb 1997). In the 
fifth instar silkworm larvae, the amount of 
SPl is higher females than in males ( Tojo et 
al. 1980; Mine et al. 1983). In the present 
study, the transcript levels of these storage 
protein genes were upregulated in both the 
infected and control larvae and reached their 
peak level at 6-9 hours after injection, alt- 
hough the transcript level in the infected 
larvae increased more sharply (Figure 1 D-F). 
This implies that a sharp increase in these 
storage proteins in the hemolymph likely con- 
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tributes to the inhibition of the generation and 
proliferation of B. bassiana at the early stage 
of infection. The increase in the transcript lev- 
el in the control larvae may contribute to 
sclerotization of the cuticle, which heals inju- 
ries to the exoskeleton caused by puncturing. 

The muscle LIM protein and chemosensory 
proteins may play important roles in the early 
events of the recovery process of skeletal 
muscles to injury and in cuticle synthesis, tis- 
sue formation, and regeneration (Nomura 
Kitabayashi et al. 1998; Stathopoulos et al. 
2002; Sabatier et al. 2003; Barash et al. 2005; 
Fort et al. 2007). LIM participation is required 
in many morphologic and histologic changes 
in insects during the transformation from larva 
to adult (Liu et al. 2007). Some chemosensory 
proteins might be involved in the immune re- 
sponse (Oduol et al. 2000; Stathopoulos et al. 
2002; Sabatier et al. 2003). In the present ex- 
periments, the transcript levels of LIM and 
CSPll were also upregulated by injecting ei- 
ther conidia or water, but the degree of 
upregulation was greater when condidia were 
injected (Figure lA and B). This result im- 
plies that muscle LIM protein and CSPll are 
not only involved in the recovery and synthe- 
sis of damaged skin but also strongly respond 
to B. bassiana invasion. However, the effect 
of the increase in the transcript levels on the 
generation of conidia and the proliferation of 
B. bassiana needs to be studied further. 

Transferrin (Tf) is a multifunctional protein 
(Nichol et al. 2002) that possesses antimicro- 
bial properties against bacteria and fungi 
(Yoshiga et al. 1997; Yun et al. 1999; Ku- 
charski and Maleszka 2003; Thompson et al. 
2003). BmTf may play an important role in 
pathogen clearance in insect innate immunity 
(Yun et al. 2009). In the present test, the ex- 
pression level of the BmTf gene increased at 
the early stage of injection and reached its 
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peak at 9 hours after injection (Figure IC). 
Perhaps BmTf affected B. bassiana gene ex- 
pression in the B. mori host cells and resulted 
in the inhibition of B. bassiana proliferation. 

Compared with our previous study on cuticle 
infection, many storage proteins and transport 
proteins were present in the hemolymph when 
the larvae were injected with B. bassiana. We 
found that the changes in the transcription 
levels of genes in the cuticle were not obvious 
in the hemolymph by real-time qPCR. These 
results suggested that the fungal pathogen re- 
sistance mechanism is different between the 
hemolymph and the cuticle. 

In conclusion, B. bassiana infection and the 
defense response in the silkworm are complex 
processes that involve signal recognition and 
transduction and the production of antifungal 
factors. Differentially expressed genes were 
identified in silkworm larvae infected with B. 
bassiana, which are probably induced by and 
responsive to the infection. These results may 
provide new insights into the complex interac- 
tions between B. bassiana and silkworms and 
may be used to gain a better understanding of 
the defense strategies elicited by the lepidop- 
teran host against fungal infection. Further 
studies are necessary to identify the functions 
of these genes and their roles in the silkworm 
defense upon fungal infection. 
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Table I. Primers used in real-time qPCR for confirmation of the differentially expressed genes. 



Gene 


Accession No. 


Forward Prmer 


Reverse Primer 


Bm GAPDH 


EU141495.1 


TTCATGCCACAAC 
TGCTACA 


AGTCAGCTTGCC 
ATTAAGAG 


arylphorin subunit alpha 
Precursor 


P14296 


TTGCCACGGTTTC 
GTTGTT 


TAGCTGCGGCTT 
CAGGATT 


chemosensory protein 1 1 


NP_00 103 7068 


GL AL 1 G 1 1 GGG 1 G 
CTGTCTAC 


TC^f^C^ A A r^f^f^f^ ATP 

1 CGGAAGGGGAl 
GTTTCTAAG 


sex-specific storage- 
protein SPl precursor 


NM_001 113276 


TCGTCGTGTAATG 
GAAAGCG 


TACCGAGGGGCA 
GCATAAGT 


transferrin 


NP_001037014 


AGTTCTAGGGCTA 
TCAGACAAAAT 


GTGACCTCTTTC 
AATGACTTCG 


muscle LIM protein 
isoform 1 


NP_00 11 03762 


CACAGCAGTCTCG 
GCATCAC 


TTATCAGCGGGT 
TTGAAAGG 


low molecular lipoprotein 
30K pBmHPC-6 (Lp-c6) 


NM_00 1044021 


AGTGCCGACTCCC 
AATACCAA 


AAGGACAAGAC 
AAGCCCGAGA 


lysozyme 


NM_00 1043983 


GCATTTCGCTGCC 
TTAGTAATGT 


TGAGGAACTGGG 
TATGTCTGGTC 


moricin 


NM_001 043364 


ACAGCCGCTCCAG 
CAAAAATACCTA 


GACCTTTACCGA 
CTGCCTTTCCTAC 
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Table 2. Differentially expressed genes from the forward SSH library and their major functional categories. 



Gene Name 


Accession 
Number 


Function Class 


Clone 
Number 


Source 


myosin 1 light chain 


AY829797 


AI Ca ion binding 


1 


Lonomia 
ohiiqua 


muscle LIM protein 
isofomi 1 


NP_00 1103762 


AI zinc ion binding 


11 


Bomhyx 
mori 


troponin C 


NPOO 1037594 


AI Ca ion binding 


1 


Bomhyx 
mori 


ferritin 


NM 001044115 


Al iron ion transport 


1 


Bomhyx 
mori 


transferrin (Tf) 


NP_001037014 


Al iron ion transport 


1 


Bomhyx 
mori 


cuticular protem RR-1 
motif 46 (CPR46) 


NMOO 1043536 


AlV cytoskeleton, 
membrane proteins 


1 


Bomhyx 
mori 


larval cuticle protein 
LCP-22 Precursor 


NP_00 1036828 


AlV cytoskeleton 
shape and membrane 
proteins 


3 


Bomhyx 
mori 


larval cuticle protein 
(LCP30 i»ene). 


X74320 


AIV cytoskeleton 
shape and proteins 


1 


Bomhyx 
mori 


beta-tubulin (Tub2) 


NM 001043499 


AIV cytoskeleton 


1 


Bomhyx 
mori 


ribosomal protein L7a 


ABS57440 


AV protein synthesis, 
ribosomal proteins 


1 


He/iconius 
mc/pomene 


ribosomal protein SI 1 
isoform 2 


NPOO 10371 11 


AV protein synthesis, 
ribosomal proteins 


1 


Bomhyx 
mori 


ribosomal protein 
S15A 


NM 001044105 


AV protein synthesis, 
ribosomal proteins 


1 


Bomhyx 
mori 


ribosomal protein SI 5 


NP 001037209 


AV protein synthesis, 
ribosomal proteins 


1 


Bomhyx 
mori 


elongation factor 1 
alpha 


ACJ71410 


AV protein synthesis 


2 


Mditaea 
casta 


translation elongation 
factor 2 (tef2), 
transcript variant 1 


NM_00 1044 128 


AV protein synthesis 


2 


Bomhyx 
mori 


H+ transporting ATP 
synthase subunit e 


NM 001098342 


AVI ATP synthesis 
enzyme 


1 


Bomhyx 
mori 


ATP synthase subunit 
c 


AFl 17583 


AVI ATP synthesis 
enzyme 


1 


Manduca 
sexta 


chaperonin subunit 6a 
zeta 


NPOO 1040 108 


AVI ATP enzyme 


1 


Bomhyx 
mori 


chitin deacetylase 


ACB54958 


AVI catabolism 
enzymes 


1 


Helicoverpa 
armi^era 


catalase 


NM 001043447 


AVII stress response 


5 


Bomhyx 
mori 


heat shock protein 1 


NM 001098297 


AVIl stress response 


1 


Bomhyx 
mori 


heat shock protein 20 


EU934239 


AVII stress response 


4 


Rpniisin 
fahaci 


alcohol dehydrodenase 


AY725189 


AVIl stress response 


4 


Arachis 


putative rubisco 
activase (Rca gene) 


AM748403 


AVII stress response 


1 


Vii^na 
nngiiicn/aia 
(cowpea ) 


imaginal disk growth 
factor 


AB 183872 


AVIIl protein 
degradation and 
processing 


1 


Bomhyx 
mori 


low molecular mass 
30 kDa lipoprotein 
19G1 precursor 


AC035751 


AIX transport protein 


9 


Bomhyx 
mori 
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Table 2. Continued. 



low molecular 
lipoprotein 3 OK 
pBmHPC-12(Lp-cl2) 


NM_001 101726 


AIX transport protein 


2 


Bomhyx 
mori 


low molecular 
lipoprotein 30K 
pBmHPC-6 (Lp-c6) 


NMOO 1044021 


AIX transport protein 


5 


Bomhyx 
mori 


low molecular 
lipoprotein 3 OK 
pBmHPC-19(Lp-cl9) 


NM 001101727 


AIX transport protein 


2 


Bomhyx 
mori 


low molecular 
iipopioiem juis. 
pBmHPC-23 (Lp-c23) 


IN IVl UU 1 1 U 1 / Zo 


AiA iranspon proiein 


Z 


Bomhyx 
mori 


peptidoglycan 
recognition protein SI 
(PGRP-Sl) 


NM 001043371 


AIX protein targeting 


1 


Bomhyx 
mori 


chemosensoiy protein 
11 


NP 001037068 


Bll intracellular signal 
transduction pathway 
molecules 


2 


Bomhyx 
mori 


lysozyme 


NM 001043983 


DI AMP 


2 


Bomhyx 
mori 


cecropin B precursor 


NPOO 103 7460 


DI AMP 


4 


Bomhyx 
mori 


moricin 


NM_00 1043364 


DI AMP 


2 


Bomhyx 
mori 


gloverin 4 


NM_00 10442 19 


DI AMP 


2 


Bomhyx 
mori 


sex-specific storage- 
protein (SPl) 
precursor 


NM_001 113276 


DII storage molecule 


3 


Bomhyx 
mori 


Arylphorin subunit 
alpha Precursor 


PI 4296 


DIl storage molecule 


12 


Manduca 
sexta 


sex-specific storage- 
protein 2 (SP2) 


NM 001044125 


DII storage molecule 


11 


Bomhyx 
mori 


putative protein 
translocase 


YP_00 1089806 


Dili metabolism and 
other proteins of 
unknown function 


2 


(lostriclium 
diftkilc 630 


27 kDa glycoprotein 
(P27K) 


NM_00 10434 13 


Dill metabolism 


1 


Bomhyx 
mori 


TPA: putative cuticle 
protein 


FAA00544 


EI; similar to known 
protein but not enough 
information to classify 


2 


Bomhyx 
mori 


putative cuticle protein 


ACY06938 


EI; similar to known 
protein but not enough 
infomiation to classify 


2 


Bomhyx 
mori 


cuticular protein 
hypothetical 43 


NPOO 1166198 


EI; similar to known 
protein but not enough 
information to classify 




Bomhyx 
mori 


unknown 




EIII; no significant 
similarity to known 
proteins 






unknown 




EIII; no significant 
similarity to known 
proteins 






unknown 




EIII; no significant 
similarity to known 
proteins 






unknown 




EIII; no significant 
similarity to known 
proteins 
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Table 3. Differentially expressed genes in the reverse SSH library and their major functional categories. 



Gene Name 


Accession 
Number 


Function Class 


Clone 
Number 


Source 


S-formylglutathione 
hydrolase, putative 


ZP_0 1750351 


AVI catabolism enzymes 


9 


Roseobacter 
sp. CCS2 


ribosomal protein 
S8 


NP_001037263 


AV protein synthesis, 
ribosomal proteins 


5 


Bombyx 
mori 


low molecular mass 
30 kDa lipoprotein 
19G1 precursor 


AC035751 


AIX transport protein 


3 


Bombyx 
mori 


putative cuticle 
protein CPH43 


ACY06938 


EI similar to known protein 
but not enough information 
to be classified 


5 


Bombyx 
mori 


unknown 




EIII with no significant 
similarity to known 
proteins 


6 
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Figure I. Variations in expression levels of fungus-responsive genes in the hemolymph of silkworms, Bombyx moh, infected with 
Beauveha bassiana. The third-day larvae of the fifth instar of Dazao strains were infected with 6. bassiana. Total RNA was extract- 
ed from the hemolymph at the indicated time points after infection and subjected to DNase I treatment and reverse transcription. 
Two microliters of each 10-fold diluted first strand cDNA (20 ng) was analyzed in each real-time qPCR reaction. The reaction was 
performed with specific primers for amplifying each of the six genes. The relative expression level of each gene at each time point 
was normalized using the Ct values obtained for the BmGAPDH amplifications run in the same plate. In each assay, the expression 
level is shown relative to the lowest expression level, which was set to one. All samples were tested in triplicate. The mean value 
± SD was used for the analysis of the relative transcript levels for each time point using the A ACt method. The 6. bassiana- 
injected and water-treated individuals are shown on the left (blue) and right (purple), respectively. A. Chemosensory protein I I; 
B. Muscle LIM protein isoform I; C. Transferrin; D. Sex-specific storage-protein SPI precursor; E. Arylphorin; F. Low molecular 
lipoprotein 30K pBmHPC-6 (Lp-c6); G. lysozyme; H. Moricin. High quality figures are available online. 
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